V 



PCT 



WORLD INTELLECTUAL PROPERTY ORGANIZATION 
International Bureau 




INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(51) International Patent Classification * : 




(11) International Publication Number: 


WO 97/18462 


G01N 23/04 


Al 










(43) International Publication Date: 


22 May 1997 (22.05.97) 



(21) International Application Number: PCT/US96/18200 

(22) International Filing Date: 25 September 1996 (25.09.96) 



(30) Priority Data: 

60/006,670 
08/591,839 



1 3 November 1995 (13.1 1 .95) US 
25 January 1996 (25.01.96) US 



(71)(72) Applicant and Inventor: WILLSON, Paul, D. [US/US]; 
89 Academy Street. Rockaway. NJ 07866 (US). 

(74) Agent: SACHS, Michael, C; U.S. Army ARDEC, AMSTA- 
AR-GCL (Building 3), Picatinny Arsenal, Dover. NJ 07806- 
5000 (US). 



(81) Designated States: AT, AU, BR. CA, CH, CN, DE, DK, ES, 
FI, GB, JP, KR, LU, MX, NO. PT. SE, SG. US, European 
patent (AT, BE, CH, DE, DK, ES, FI, FR, GB. GR, IE, IT, 
LU, MC, NL, PT, SE). 



Published 

With international search report. 
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(57) Abstract 

An apparatus and method for automatic recognition of concealed objects and features thereof, such as contraband in baggage or 
defects in articles of manufacture, is disclosed. The apparatus uses multiple energy x-ray scanning to identify targets with a spectral 
response corresponding to a known response of targets of interest. Detection sensitivity for both automatic detection and manual inspection 
are improved through the multiple-energy, multispectral technique. Multichannel processing is used to achieve high throughput capability. 
Target identification may be verified through further analysis of such attributes as shape, texture, and context of the scan data. The apparatus 
may use a statistical analysis to predict the confidence level of a particular target identification. A radiograph, CT image, or both may be 
reconstructed and displayed on a computer monitor for visual analysis by an operator. Finally, the apparatus may receive and store input 
from the operator for use in subsequent target identification. 
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APPARATUS AND METHOD FOR AUTOMATIC RECOGNITION OF 
CONCEALED OBJECTS USING MULTIPLE ENERGY COMPUTED TOMOGRAPHY 

GOVERNMENT INTEREST 
The invention described herein may be manufactured, used and licensed by or for the 
government for governmental purposes. 

BRIEF DESCRIPTION OF THE INVENTION 
The present invention relates generally to the field of nondestructive inspection. More 
particularly, the present invention relates to an apparatus and method for detecting concealed 
objects and features thereof, such as contraband in baggage or defects in articles of manufacture, 
using multiple energy computed tomography. 

BACKGROUND OF THE INVENTION 

Conventional x-ray scsmning is used in a number of fields to detect objects or features not 
visible to the human eye. For example, in the medical and dental fields, x-ray systems are used to 
detect features of interest in rendering a clinical diagnosis, such as a fractured bone or a cavity. In 
the manufacturing industry, x-ray systems are used similarly to inspect parts for defects 
Fractures or voids below the surface of a weld, for example, can be detected from an x-ray 
image, thus avoiding possible failure of the part should it be used in its defective condition X- 
ray systems are also used in airports and other public facilities to inspect containers for weapons, 
explosives, and other contraband. 

In each of the foregoing applications, the x-ray system is merely an imaging device 
without the capability of automatic identification of targets. These systems produce a gray scale 
image, representing the total x-ray energy absorbed by all objects between the x-ray source and 
the image plane: the more energy absorbed, the darker the corresponding spot on the image 
Usmg this projection method, the resulting images or radiographs are often difficult to interpret 
because objects are superimposed. Data obtained from x-ray images are generally unsuitable for 
automauc detection because of the complexity involved in resolving superimposed objects A 
trained operator must carefully study and interpret each image to render an opinion on whether or 
not a target of interest is present Where an application requires a large number of radiographs to 
be interpreted, operator fatigue and distraction can compromise detection capability 

X-ray Computed Tomography (CT) is a technique that produces an image of a cross- 
sectional slice of an object from a series of attenuation measurements taken at various angles 
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around the object. The CT image does not suffer from the superpositioning problem presented 
with standard radiographs. Although CT data can provide precise, quantitative information about 
the characteristics of objects in the scan plane suitable for automatic detection of targets, it too 
has limitations. Conventional CT systems take considerable time to perform a scan, to capture 
the data and reconstruct an image. Throughput of CT systems is low. Coupled with the size and 
expense of conventional CT systems, this limitation has hindered CT use in applications such as 
baggage or parts inspection where object throughput is a major concern. 

U.S. Patent 5,367,552 to Peschmann (1994) teaches one method for improved CT 
throughput. In the Peschmann system a conventional x-ray scanner is first used to pre-scan an 
object, followed by CT scanning at locations selected from analysis of the pre-scan data. 
Although the solution taught by Peschmann provides improved detection capability over 
conventional x-ray systems, it has several limitations. First, it requires pre-scanning of the object 
with a conventional x-ray system which takes time and provides limited results as discussed 
above. Second, in order to save time, a CT scan is performed only at selected locations which 
could result in failure to identify targets of interest, especially where the target is masked or 
otherwise difficult to detect with a conventional x-ray scanner. Third, because the Peschmann 
invention uses a conventional rotating CT device, the throughput is limited by the mechanics of 
the rotation. Fourth, the baggage is stopped, again limiting throughput, to allow for rotation of 
the x-ray source around the object at each slice. Finally, Peschmann teaches the use of 
conventional single- and dual-energy techniques for generating CT data whereas a multiple- 
energy or multispectral technique as described herein would result in improved target 
identification. 

OBJECTS AND SUMMARY OF THE INVENTION 
25 Accordingly, the objects of the present invention are as follows: 

(a) to provide automatic recognition of concealed objects and features thereof with or 
without operator involvement; 

(b) to detect small quantities or features of a target of interest concealed within an object; 

(c) to aliow for high throughput of objects during scanning operations without 
30 compromising detection capability; 

(d) to provide CT data using a compact stationary x-ray source array and detector array; 

(e) to provide an enhanced x-ray image or CT image or both for operator viewing; 

(f) to provide statistically based confidence levels for target identification; and 
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These and other objectives are achieved by an apparatus and meUtod for aulomatic 
recog„,„o„ of concealed obje«s. The apparatus uses multiple energy x-ray uJZ^ .•«. 
targets with a spectral response corresponding to a known Ksponse^f^argeKofinterest ^ 
Detectton sensitivity for both automatic detection and manual LyauZT- ^ 
^.-gy.multis^^ 

^nghpm capabUhy. Targe, identification may be verified through Lher 

aanbures as shape, texntre, and context of the scan data. The apparams may use a lIZl 

■mage, or bom may be reconstructed and displayed on a comnuter moni„> i Tf* CT 
an operator. Finally, tire apparahts may receive Ld ZlZTon Z " T "** * 
subsequent urge, identification. ^ *' ° Pen " 0r for in 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram of the individual hardware components of the ore»n, 
invem.on, with their logical connections. ts ot the present 

" <"~ M,b ° di " K " 1 ° f *< Mention for 

FIG. 4 is a schematic representation of another embodiment of the present invent r 
inspection of manufactured parts. P invention for 

HG 2 " 3 \ hart dCPiCting "* maJ ° r StePS i0 out the present invention 

FIG. 6 is a flow chart depicting the major processing steps on the present i^Z r 
analyzing multispectral data. H> «n me present invention for 

^^.Tis an illustration <^*v»«n*^^^ na ^ 

DETAILED DESCRIPTION OF THE INVENTION 
Referring to FIG. 1 , an apparatus in accordance with th~ „~ 
»ay of overview, me apparahrs Is firs, descrihed^Z^et EZ, ™ m ""' " **** * 
parricular embodiment source array 1 is an L shin ' <°n>ponen B . In tins 
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FIG. 1 , provide a series of co-planar fan beams 16-1 9 in response to a signal from controller 6. 
Each fan beam 16-19 is comprised of x-ray photons of varying energy levels within a fixed 
energy band. Source amy 1 is positioned opposite to and separated from detector array 2 to form 
a space within which object 21 may be scanned. 

Detector array 2 is also L-shaped in this particular embodiment and is comprised of 
individual detector elements 20, which are arranged at regularly spaced intervals along its length. 
Detector element 20 absorbs the photons from source array 1 and provides a voltage signal to 
data acquisition circuit 3. The resulting voltage signal is proportional to the energy level of the 
photon absorbed by detection element 20 after traveling along the beam path from source airay 1 
to detector array 2. By using a series of separate detector elements 20, the position along detector 
array 2 where the photon was absorbed can be determined and correlated to the corresponding 
energy measurement. The spatial resolution of the resulting x-ray scan data is thus determined by 
the size and spacing of detector elements 20. Very small detector elements 20 closely arranged 
and abutted against one another along array 2 will provide precise information as to where a 
photon energy measurement was taken. Such a configuration, however, requires a large quantity 
of elements 20 to fill the length of array 2 and a correspondingly large number of data acquisition 
circuits 3, as described below, to capture and process the resulting data. Thus, an analysis of the 
required spatial resolution versus component cost should be made to determine the optimal size 
and spacing of detector elements 20 along array 2. 

A multichannel data acquisition circuit 3 is provided for rapid processing of the voltage 
signals from detector elements 20. Each channel of circuit 3 is electrically coupled to one 
detector element 20 and is comprised of a detection processor 1 1 and a series of comparators 9 
and counters 10. Detection processor 1 1 provides a separate voltage threshold signal to each 
comparator 30-34 for comparison with the voltage signal from detector element 20. As more 
fully described below, data acquisition circuit 3 forms spectral attenuation data, defined as the 
number of photons absorbed by each detector element 20 versus the energy level of the absorbed 
photons, from the analog voltage signals provided by detection element 20. These spectral 
attenuation data are provided to processor 4 at specified intervals for further processing. 

Processor 4 is a computer processor with parallel processing capability. Each channel of 
data acquisition circuit 3 is electrically coupled to processor 4 for input of spectral attenuation 
data. Processor 4, using CT reconstruction algorithms, reconstructs from the spectral attenuation 
data a tomographic image of the cross-sectional slice of object 20 which was scanned to acquire 
the data. A pixel in a tomographic image of a cross-sectional slice is called a voxel because it 
represents material in a volume element whose dimensions are those of the pixel and the 
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titickness of the die. Tie raw x-ray image dara are als. preserved by processor 4 for s „. 

each of the five measured attenuation data seta resulting in multiple-cnerg, CT dam As L, 

data. He resulting CT .mage dara rue then matched with the data from file server 5 which 
represent kno w, targets of interest. Objects or features of interest tha, m c^ ^ 
object space 2. are thus automatically identified. Selected romographic images and raw data are 
also stored on file server 5 for later reference. ruwoaraare 
Tagging system 8 is provided in this particular embodiment for applying an identification 
tag toobjectj. infuse toa signal fiom con*o ller 6 who, a p^tive^X^oT 
made by the appara.ua. A defective pa* or baggage tha, contains commons Z 5l 
^-<»™en.ion.na^^^^ 

handhng system may be used ,„ automatically separate defective parfs from noLective p^ 
U -'"^*P™tdeJ t odi S p.ayenh^^^^ 

Z^lfT**" *" * h<mm °<™°'- 7 iselecuicaHy coupLT^ 

ZT T"Z "f ^ *"• 7 ^ iS — '° seLed images 

from file server 5 for display to the human operator X-rav and rr i m ™ .. ™' ma 8es 

-*<<•*• *• 'tat .n the preferred embodiment, interface 7 is a^ge.Wgh 
resolution color touchscreen for ease of viewing and input operation 

Having provided a general overview of the major components, attention is now turned to 
a mora detaded operational description. Depending upon the particular appUcZ"™ 
belt or other means is provided to move object 2, forward through die space defin d by 3 
amy I and detector anay 2 for scanning. As object 2 1 is moved fonva^ each of Z x-ray 
sources 2-,5 isactivated in seouence by controller 6 such that only one x-ray J „ 
produced a. any one inaant. The path of a particular photon lha, caused a volLe sienal K 2 

„h„,„ ^ 0taampk - *°* S0UTO 12 is »i™«< f.«t, emitting a continuum of enerny 
Ph^ns of ^ spectral conlent „ „ Mme ^ enKgy Uvei 

wfcehpassesthroughsite 22 in object 21. As each photon is absorbed by debtor eZen't^a 
vo -age signa, in proportion ,o the energy level of tire photon is provided ,„ c,r , ' 
tar Cement. The voltage signa! is provided as a first input signa, to each of! five 
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comparators 30-34. Each comparator 30-34 receives a second input signal from detection 
controller 1 1 for use as a threshold voltage. When the input voltage from element 20 exceeds the 
threshold voltage from controller 1 1 , a counter 35-39 is incremented by one unit For example, 
when the input voltage of comparator 30 exceeds the threshold voltage from controller 1 1, 
5 counter 35 is incremented by one unit Comparators 30-34 and counters 35-39 function in the 
same manner to count the number of times a threshold voltage is exceeded. In this manner the 
total count in counter 35 accumulated while source 12 is active represents the intensity of 
spectral range above the comparator threshold. Intensity of multiple spectral ranges are thus 
formed by setting the voltage thresholds of comparators 30-34 at different levels. After a brief 

10 period of time, x-ray source 1 2 is deactivated by controller 6, the contents of counters 35-39 are 
loaded into processor 4, and the counters are reset to zero in response to a signal from controller 
6 via detection controller 1 1 . X-ray source 1 3 is then activated by controller 6 and the detection, 
data acquisition, and processing steps described above are repeated using source 13, followed by 
source 1 4, and finally, source 15. Double buffering is employed between counters and controller 

15 6 for greater throughput. The counter data are sent to processor 4 for the period during which 
each of the x-ray sources was activated is sufficient to calculate one set of tomographic images. 

The use of a parallel data acquisition circuit 3 for each detector element 20, increases the 
processing speed of the apparatus so as to limit the amount of time needed to activate source 12 
and collect the resulting attenuation data. In this way, object 21 may be moved continuously 

20 through the apparatus during scanning operations while still providing sufficient data for 
automatic target identification. Parallel data acquisition is made practical through use of 
Application Specific Integrated Chips (ASIC) for circuit 3. Data acquisition circuit 3 also 
provides for field adjustment of the selected energy ranges. 

An alternative energy resolving method would be to send the photon signal to an analog 

25 to digital converter (ADC) followed by digital binning. For a limited number of bins the first 
embodiment is faster, cheaper, and preferred because digital processing is avoided and energy 
discrimination is done in the analog mode. If the number of bins were increased to several 
hundred, then multiple ADCs would be more practical. 

As fan beam 1 6 passes through object 2 1 , some of the photons of beam 1 6 are absorbed 

30 by the material at site 22 within object space 21, some photons pass through unaffected, and 
some photons are scattered showing up as lower energy photons. In that the path of scattered 
photons is unknown, they are a source of noise. A collimator placed in front of detector element 
20 decreases the portion of scattered photons impinging the detector element. The number of 
photons absorbed compared to the number of unaffected photons, that is the attenuation of the x- 



6 



10 



15 



20 



25 



30 



WO 97/18462 PCT/US96/18200 
rays, is a function of the energy of the photons and the type of material at site 22 through which 

20^1 n 8 T * DUmber ° f Ph ° t0nS ° f ^ Cnergy 13086 « ^tector element 

20 will be less after passing through the material of object 21 than would have been if no 

material were present. In as much as the proportional decrease in number of photons of different 

energy ,s a function of the chemical composition of the material, the ratios representing the 

proportional decrease can be used to characterize the material. By using sources 12-15 having a 

contmuum of energy, and capturing the resulting multi-energy spectrum after the photons have 

would be captured with no materia, present atsite22, 

chemical composition of object 2 1 over systems that use a single- or dual-energy technique 
Used in conjunction with other image processing means, such as described below this * 
multispectral data provides for greater discernment of contraband in baggage defects in a «i e l« 
of manufacture, and similar applications. 8 m 

Processor 4 is used to separate the counter values by spectral range and to calculate a 
tomograph* image for each spectral range. Attenuation of the x-rays by the material at site 22 is 
contained in the data acquired from those detector elements which acquired x-ra s p^ng 
through site 22. Tomographic reconstruction of the voxel at site 22 is performed by appropriate 
mathemancalcombinationofthedata ^B^cctB^^^^^ 
x-ray detector of figure , the algebraic reconstmction technique is most appropriate for 

calculating the tomographic image. For other geometries other reconstruction techniques may be 
n^ 

20 at the same time, the voxels in the respective spectral* differing tomographic images are abo 
"Thesetofvaluesformefive 

the actual spectral attenuation occurring at that voxel and can be used for characterizing the 
material at the voxel site. For object 2 1 in continuous motion while attenuation data are acquired 
bydetector elements 20, the scan ishelical and the tomographic reconstruction is a type o^ 
hehca tomographic reconstruction. Processor 4 after calculating the vahies stores upvalues 
onto file server 5 for later reference. 

P ""^ r4fiKte ^^^»«<»feachvoxel.„taofvario U skno TO ,a,Betvalue S 
stored on fi,e server 5. From the conspanson e coded itnag, is e„„ smK .ed which exem^T 

tal LIT '" *•* «" "*"* * «*— • » - set of 

™ag.s from the d.fference tn the values of the respective voxels in the respective spectral 

-ages. ,„ the, case the difference images are f„ ,„ urge, values instead cL stages 
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themselves to construct a coded image which exemplifies voxel-by-voxel the most probable 
materials and the degree of probability. Processor 4 performs an analysis on the coded images 
by looking for groupings of neighboring voxels whose composition is nearly the same. Processor 
4 compares the grouping size and shape to that of the probable material to further uniquely 

5 identify the object Processor 4 may be used to perform other analyses to measure or calculate 
unique parameters which are known to increase the probability of recognizing objects. The data 
can be Fourier transformed and a spatial frequency characterization made. The data, for 
example, can be wavelet transformed to enhance finding object shapes. Characteristics such as 
3-dimensional roundness, granularity or texture can be measured. These and other data 

0 processing techniques are well known in the art and may be practiced with the present invention 
for certain applications. The final result of this data analysis can be either a decision on whether a 
particular feature was detected, or can be an enhanced image for human viewing, or both. 

Because x-ray systems can drift considerably from the desired calibration state, the 
present invention also includes means for reducing this problem. A gross calibration factor can 

5 be added in real-time to the overall system using processor 6 along with a standard in the field of 
view. In this way, a general calibration correction is applied to all of the detector elements 20. 
The calibration drift problem is also reduced in the present invention through the use of relative 
absorption values, which serve to increase the signal-to-noise ratio of the data as described 
below. 

0 Tagging system 8 is used in one embodiment to automatically tag a container or part 

when a positive identification is made. An alternative embodiment could include an automatic 
sorting and material handling system in place of tagging system 8. User interface 7 is used to 
both display the results for human viewing and to receive input from the operator, initial 
calibration of the system can be achieved using interface 7 to input the correct response for scan 

5 objects of known content. The same feature can be used for continuous improvement of the 
system. 

Continuous improvement is performed by a self-learning method provided for in the 
system. Objects 21 of known material and identification are scanned by the system. The system 
in performing its feature discerning operations will find voxels of like material, group the voxels 
o into shapes, calculate parameters characterizing shape, size, and texture even though the new 
object does not correlate to any known object. The operator groups the calculated parameters 
and tells the system to add the new set of parameters to its list of known objects. The system 
now includes the new object in its matching for automatic detection. 
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Having provided a general overview of the present invention, attention is turned to a 
more deta.led description of suitable physical components for certain applications FIG 2 is a 
schematic representation of one embodiment of the invention as developed to meet the specific 
needs of a particular baggage inspection application. In this application, baggage throughput as 
well as system cost and size axe major concerns. The system must function at normal baggage 
handhng rates (up to 2 feet per second), fit in a limited space, and must be compatible w*n 
existrng baggage handling equipment. This is accomplished with the apparatus illustrated in FIG 
2. An array of x-ray sources 40 and photon energy absorbing detectors 41 are arranged in a 
verucal plane between the ends of two conveyor belts such that a bag moving from one belt to 
the other w,ll come between the x-ray source array 40 and detector array 4 1 The spacin fi of 
sources and detectors provide multiple beam paths, which allow for computed tomography with 
no movmg parts within the apparatus. Conventional CT systems move the x-ray source the 
detectorarmy.orbothtoachievesimilarresults. Although other arrangements can be used the 
square geometry of the source and detector arrays provides sufficient room for baggage to pass 
through, while minimizing space requirements for the apparatus. Linear and relatively short 
components allow foramodular design ^ can be adapted to various geometries. Tbe arrays 40 
and 4! are located at a gap in the conveyor system to reduce signal interference from the 
conveyor hardware. Finally, the system may be designed such that it can be scaled geometrically 
in stze and resolution while still employing the same mathematical analysis and using the same 
baste components. In this way, systems can be easily customized for a particular application 
without a huge design effort. ppucauon 

FIG. 3 is a schematic representation of another embodiment of the source and detector 
■ekments suitable for an application where the size of the apparatus is not a major concern. In this 
embodiment, the source array 52 and detector array 50 are configured as two concentric rin s of 
equal d.ameter with the counterpoint of one array offset longitudinally a short distance from the 
counterpomt of the other. The distance between source array 52 and detector array 50 is 

d vantage of broader coverage for each source in comparison with the first embodiment of FIG 
because of the larger distance between source and object Also, the scan space is evTnly 
covered by x-ray beams in this embodiment because of its symmetry. The major disadvantage of 
*» embod.ment is that it occupies more space than the L-shaped arrays of FIG. 2. It JZL 

wuh a senes of hnear arrays to form a polygonal shape. In this way, a modular design can be 
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achieved that approximated the circular shape of FIG. 3 while reducing component replacement 
costs. 

In other applications, particularly the manufacturing arts, the object to be scanned may be 
uniform in shape, composition, or both. True tomography may be replaced with pseudo 
5 tomography by employing a priori knowledge of manufactured components-their geometry, 
position, and material composition-as constraints in the reconstruction. In such situations, 
reduced system cost and faster computational speed can be achieved using a source and detector 
embodiment as illustrated in FIG. 4. Here, two sources 60 and 61 are used with two detector 
planes 62 and 63. Although greatly simplified in comparison to the other embodiments of FIG. 1- 

10 3, the techniques of the present invention can be used in this and similar embodiments. 

Having provided an overview of suitable hardware elements to be used in accordance 
with the present invention, attention presently turns to the processing steps to be executed on the 
hardware. Turning to FIG. 5, the major steps in carrying out the present invention are illustrated 
as a flowchart. The process begins with initiation of the source number to one in step 80. The 

1 5 source is activated in step 8 1 to generate a beam of photons directed towards the object to be 

scanned. The energy level of the photons after traveling along the beam path is then measured in 
step 82. In step 83 spectral data are formed from the measured energy levels by counting the 
number of photons within a certain energy range absorbed by each detector during a specific time 
increment. In step 84 the source number is checked to determine if all sources have been 

20 activated. If not, process step 88 is performed to increment the source number and return to step 
81 . If all sources have been activated, then process step 85 is performed to analyze the spectral 
data. In this way, steps 81-83 are repeated for each x-ray source. 

The results from the above steps are used to form and analyze multispectrai CT data in 
step 85 as discussed fully below. These data are then compared to data representing known 

25 targets of interest in step 86 to determine if one or more targets of interest are present. If not, an 
appropriate message is displayed in step 89 along with an enhanced x-ray or CT image for 
operator viewing. If one or more targets are identified in step 86, then the container is tagged or 
sorted in step 87 and an enhanced image is displayed in step 90 with a unique color or texture to 
identify the targets along with an appropriate text message to the operator. 

30 The processing steps used to identify targets of interest are illustrated in greater detail in 

FIG. 6. A CT algorithm is applied to each of the spectral data sets in steps 101, resulting in 
multispectrai attenuation data for each voxel comprising the object space. Calibration corrections 
are applied in step 102 to correct for the variation among detector elements and the calibration 
drift that occurs over time. This step may be performed before CT reconstruction, after, or both 
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depending on the level of real-time calibration requtted for .he particular application. Use of a 
standard m ore Held of view geKrally requires the ealibration eorrection be applied after CT 
reeonstmction in order .0 detennine the p K) per coneclion values. The fixed variation among 
de.ec.or element ^suiting from imperfections in one manufacturing process on tt,e other hand 
cm be coveted before CT reconsttaction by applying a correction direcUy ,o fte signal outou. 
of each detector element. v 

After CT reeonstntolion and cal.bra.ion correction, one spectral image is combined witii 
anoti^r u, step .03 to obtain mlative attenuation date. This operation improves me signal-to- 
no.se ra..„ of me resulting date. A malched filter is applied to me resulting dato se.s in step 104 
.0 determine me stetistical likelihood orprobabiltty of a ma«ch .0 tergos of taction a voxel- 
by-voxe bas,s. These dam may be rep.esen.ed by me probabilhy ma. a particular voxe. comains 
a part.cular material detemtined b, comparing me measured aquation dati, with attenuation 
date prevoualy measurod for tergete of interast Because mis comparison is made for «*h of 
several targete each voxe. may have severe, possible matehes witi, va*ng .evels of probability 

In step 105 contiguous or near contiguous voxels with similar matches are linked for 
toher analys* beyond me voxel-by-voxel comparison of step 104. Shape, size, textore, and 
ofter character.^ of me date obtoined from step 1 05 are matehed .0 similar date reprasenting 
targets of .uteres, instep ,06 to provide grea.er.arge. discertunen,. ,„ tin, way, me £ZZ 
level of a n^ch , s lnOTaMd ^ to ^ of for example ^ 

qualiues of.be items to be identified. 

I " s ^>08te"lrtveattenuationdaufrt ) ms.epl03isfusedorcombinedtof„rma 
angle ,mage ,or operator viewing. This image is enhanced in step 1 09 witi, infonuation ftom me 
ansjys.s of steps 104-106. For examp,e, highly eortelated voxels can be illusttated in JZ g T 
wtth a common color or texmre to assis. in operator viewing of toe intonation. The result of 

~ t"^ " 3 ^ fa "» ' 07 M M 1 ^PWc image 

m step 1 10, wh.ch may be performed siroulteneously on the same user interface screen 

Turnrog to HO. 7. me processing of multispecteal date to fonn an enhanced image whh a 
greater s^al-to-noise ratio is illusttated. images 71-74 corrospond .0 tire five energy ra»geT 

ttZTjd 5 ^ 7 °" 74 *™ 0bjK,S "~ * > remand 
70 7* K k" T rePrcSa,ttd " y 3 *» H,K a "«- in all five images 

*W4. tire otoer tome objecte appear wift varying intensities. This simulates a key differeT 

etZTr ^ ^ — " * - «* -i- 0^.0, 

— , ^ defM ' iVe de,mM e ' Cmra,S ' md of .he 

apparams and processus metimds. .„ many cases, tire attenuation values of acmal object 
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the scan space will vary with intensity level of the x-ray source while the artifact will have an 
attenuation value independent of energy level. The circle, rectangle, and triangle shown in 
images 70-74 simulate the variation of attenuation that different materials can have with varying 
x-ray source intensities. Thus, the circle is best resolved at the energy level of image 70, while 
the rectangle is best resolved at the energy corresponding to image 72 and the triangle at the 
energy level of image 74. 

The fusing process in the present invention is used to both eliminate the artifact 
represented by the line in images 70-74 and to produce a single enhanced image 75 that includes 
a clear picture of all three objects. The first step in this fusing process is to subtract one of the 
images, say 71 , from the other four images 72-74 to eliminate artifacts in the image data. 
Because the artifacts, as represented here by the line, do not vary in intensity from one image to 
another, this differencing technique effectively eliminates the artifacts from the resulting images. 
These four images are then summed to produce an enhanced image 75 that includes unique 
colors, numerical values, or other distinguishing qualities assigned to identify the objects. Since 
the multispectral data are collected electronically by the data acquisition circuit 3 of FIG. 1, this 
process is performed numerically in the present invention. The final image 75 may be used to 
display results to an operator and for subsequent processing, including shape, wavelet, fractal or 
other techniques of image data analysis. 

A number of advantages of the present invention are evident from the above description. 
First, the invention provides a means and apparatus for automatic detection of concealed objects 
with or without operator involvement. Small quantities or features of a target of interest 
concealed within an object may be detected. The invention provides for high throughput of 
objects during scanning operations without compromising detection capability. CT data are 
obtained using a compact and stationary x-ray source and detector array. An enhanced x-ray 
image, CT image, or both are provided for operator viewing. Statistically based confidence levels 
for target identification may be used based upon the data stored within the system, and a 
continuous learning capability is provided for improving target identification with system use. 

The foregoing description of specific embodiments of the present invention have been 
presented for purposes of illustration and description. They are not intended to be exhaustive or 
to limit the invention to the precise forms disclosed, and obviously many modifications and 
variations are possible in light of the above teachings. For example, the number and arrangement 
of x-ray sources and detectors can vary considerably depending on the application. Where space 
constraints are not a factor, the sources can be moved away from the detectors resulting in greater 
coverage per source thus reducing the number of sources required. Also, the number and 
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sophistication of da, pressing Heps cm va* greatly depending on the ^ identification 
resolunon needed fora panicularappfication. CT processing, for examp,.. may be t^neeeasary 
where objects ,n the sear, plane are relatively thin and homogeneous. The same holds tnre fi7 
shape, stze and texture analysis: in some applications a voxel-by-voxel comparison wil. provide 
sufficten. dtscemmen. capability, in omer applications shape, size, or rexttue analysis may be 
required to meet target identification requirements. 

The embodiments illustrated and described above were thus chosen to best explain the 
pnncples of the mvenrton and its practical application, thereby enabling outers skiHed in the an 
" be* unto the mventton and various embodiment with various modification as am sulTo 
0, pamctdar use contemplaKd. „ is int.rn.ed ma, the scope of the invention be defined by tit 
following Claims and their equivalents. y 



13 



WO 97/18462 PCTAJS96/I8200 
I claim: 

1 . An apparatus for detecting targets of interest concealed within an object, said apparatus 
comprising: 

an x-ray source for producing a beam of photons having an energy spectrum within a set 
5 band, said beam directed towards said object; 

a photon energy absorbing detector array positioned opposite said x-ray source such that 
said object is located between said source and said detector anay for measuring the energy 
attenuation of said photons after traveling along said beam path and through said object; 

a data acquisition circuit responsive to said detector for forming attenuation spectral data 
l o from said photon energy measurements; and 

means for analyzing said attenuation spectral data to identify targets of interest within 
said object 

2. An apparatus as recited in claim 1 , wherein said data acquisition circuit includes a 
plurality of comparators and counters. 

15 3 - An apparatus as recited in claim 2, wherein said analysis means includes computed 

tomography. 

4. An apparatus as recited in claim 3, wherein said analysis means further includes the 
application of calibration coefficients. 

5. An apparatus as recited in claim 3, wherein said analysis means further includes a 
20 mathematical operation to form relative attenuation data. 

6. An apparatus as recited in claim 3, wherein said analysis means further includes 
application of a match filter to compare measured attenuation values with those of said targets of 
interest. 

7. An apparatus as recited in claim 5, including a means for fusing said relative 
25 attenuation data for display. 

8. An apparatus as recited in claim 6, wherein said analysis means further includes the 
linking and analysis of voxels having similar attenuation values. 

9. An apparatus as recited in claim 2, wherein said x-ray source is an L-shaped array 
comprised of a plurality of individual x-ray sources. 

30 10. An apparatus as recited in claim 2, wherein said detector is an L-shaped array 

comprised of a plurality of individual x-ray detector elements. 

1 1 . An apparatus as recited in claim 2, wherein said x-ray source is a circular array 
comprised of a plurality of individual x-ray sources. 
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1 2. An apparatus as recited in claim 2, wherein said detector is a circular array comprised 
of a plurality of individual x-ray detector elements. 

13. An apparatus as recited in claim 2, wherein said x-ray source is composed of a set of 
linear arrays, said linear arrays comprised of a plurality of individual x-ray sources. 

5 14. An apparatus as recited in claim 2, wherein said detector is composed of a set of 

linear arrays, said linear arrays comprised of a plurality of individual x-ray detector elements. 

1 5. An apparatus as recited in claim 8, including a means for teaching the system a new 
set of parameters to be used in subsequent target identification. 

16. A method of detecting targets of interest concealed within an object, said method 
10 comprising the steps of: 

generating a beam of photons having an energy spectrum within a set band, said beam 
directed towards said object; 

measuring the energy level of said photons after traveling along said beam path and 
through said object; 

forming attenuation spectral data from said photon energy measurements; and 
analyzing said attenuation spectral data to identify target data within said object 

17. Tbe method of claim 16, wherein said analysis step includes computed tomography. 

1 8. The method of claim 16, wherein said analysis step includes the application of 
• calibration coefficients. 

1 9. The method of claim 16, wherein said analysis step includes mathematically 
combining one attenuation data set with another attenuation data set to form relative attenuation 
data. 

20. The method of claim 16, wherein said analysis step includes the application of a 
matched filter. 

2 1 . The method of claim 20, wherein said analysis step further includes the linking and 
analysis of voxels having similar attenuation data. 

22. The method of claim 19, further comprising the step effusing said relative 
attenuation data to create an enhanced image for display. 

23. The method of claim 22, wherein said enhanced image includes a unique display of 
30 identified targets to aid operator viewing. 

24. Tbe method of claim 21, further comprising the step of teaching the system a new set 
of parameters to be used by the system to identify new objects. 
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